o 
o 

(N 

C 
G 






Luttinger liquid versus charge density wave behaviour in the 
one-dimensional spinless fermion Holstein model 

H. Fehske a '*, G. Wellein b , G. Hager b , A. Weifie c , K. W. Becker d , A.R. Bishop 

a Institut fur Physik, Universitdt Greifswald, D-174-S7 Greifswald, Germany 

Regionales Rechenzentrum Erlangen, Universitdt Erlangen-Niirnberg, Germany 

c School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia 

Institut fur Theoretische Physik, TU Dresden, D-01062 Dresden, Germany 

6 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 



T3 

C 

o 
o 



> 

CO 
CN 

O 
^D 
O 
^1" 

O 

-i— > 

a 



T3 

C 
O 

o 



x 



Abstract 

We discuss the nature of the different ground states of the half- filled Holstein model of spinless fermions in ID. In the 
metallic regime we determine the renormalised effective coupling constant and the velocity of the charge excitations 
by a density-matrix renormalisation group (DMRG) finite-size scaling approach. At low (high) phonon frequencies 
the Luttinger liquid is characterised by an attractive (repulsive) effective interaction. In the charge- density wave 
Peierls-distorted state the charge structure factor scales to a finite value indicating long-range order. 
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The challenge of understanding quantum phase tran- 
sitions in novel quasi- ID materials has stimulated in- 
tense work on microscopic models of interacting elec- 
trons and phonons such as the Holstein model of spin- 
less fermions (HMSF) 
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The HMSF describes tight-binding band electrons 
coupled locally to harmonic dispersionless optical 
phonons, where t, u>o, and g denote the electronic 
transfer amplitude, the phonon frequency, and the 
electron-phonon (EP) coupling constant, respectively. 
Despite of its simplicity the HMSF is not exactly 
solvable and a wide range of numerical methods has 
been applied in the past to map out the ground-state 
phase diagram in the g-wo-plane, in particular for the 
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half-filled band case (N e i = N/2). There, the model 
most likely exhibits a transition from a Luttinger liquid 
(LL) to a charge-density wave (CDW) ground state 
above a critical EP coupling strength g c {wo) > [1]. 

In this contribution we present large-scale DMRG 
calculations, providing unbiased results for the (non- 
universal) LL parameters u p , K p , and the staggered 
charge structure factor S c (tt). 

To leading order, the charge velocity u p and the cor- 
relation exponent K p might be obtained from a finite- 
size scaling of the of the ground-state energy of a finite 
system Eo(N) with N sites 

eo(oo) - (E (N)/N) = (ir/3)(u p /2)/N 2 (2) 

(eo(oo) denotes the bulk ground-state energy density) 
and the charge excitation gap 

A ch (iV) = E^i^-EoiN) = tv(u p /2)/(NK p ) (3) 

(here Eq(N) is the ground-state energy with ±1 
fermions away from half-filling) . The LL scaling rela- 
tions (2) and (3) were derived for the pure electronic 
spinless fermion model only [2] . 
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Fig. 1. Finite-size scaling of the charge gap A C ^(AT) and the 
ground-state energy Eo(N). Exact diagonalisation (ED) data 
is included for comparison. 

Figure 1 demonstrates, exemplarily for the adiabatic 
regime, that they also hold for the case that a finite EP 
coupling is included. The resulting LL parameters are 
given in Tab. 1. Most notably, around uio/t ~ 1, the LL 
phase splits in two different regimes: For small phonon 
frequencies the effective fermion-fermion interaction is 
attractive, while it is repulsive for large frequencies. In 
the latter region the kinetic energy is strongly reduced 
and the charge carriers behave like (small) polarons. In 
between, there is a transition line K p = I, where the 
LL is made up of (almost) non-interacting particles. 
The LL scaling breaks down just at g c (ui/t), i.e. at 
the transition to the CDW state. We found g'^iuj/t — 
0.1) ~ 7.84 and gl{io/t = 10) ~ 4.41 [3]. 
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Table 1 

LL parameters at small and large phonon frequencies. 

Figure 2 proves the existence of the long-range or- 
dered CDW phase above g c - Here the charge structure 
factor 
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unambiguously scales to a finite value in the thermo- 
dynamic limit (N — > oo). Simultaneously A c h(oo) ac- 
quires a finite value. In contrast we have S c (tv) — > in 
the metallic regime (g < g c ). The CDW for strong EP 
coupling is connected to a Peierls distortion of the lat- 
tice, and can be classified as traditional band insula- 
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Fig. 2. Scaling of the charge structure factor S c (tt) using pe- 
riodic (PBC) and open (OBC) boundary conditions. 

tor and bipolaronic insulator in the strong-EP coupling 
adiabatic and anti-adiabatic regimes, respectively. 

The emerging physical picture can be summarised 
in the schematic phase diagram shown in Fig. 3. 
In the adiabatic limit (u>o — * 0) any finite EP cou- 
pling causes a Peierls distortion. In the anti-adiabatic 
strong EP coupling limit {luq —* oo), the HMSF can 
be mapped perturbatively onto the XXZ model and 
the metal-insulator transition is consistent with a 
Kosterlitz-Thouless transition at gl{oo) ~ 4.88. 
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Fig. 3. Phase diagram of the ID half-filled spinlcss fcrmion 
Holstcin model. Here a — cJo/i, 6 P — uJog , and A — e p /2t. 
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